Searching a Database under Decoherence 
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We report on the effects of a simple decoherence model on the quantum search algorithm. 
Despite its simplicity, the decoherence model is an instructive model that can genuinely imitate 
realistic noisy environment in many situations. As one would expect, as the size of the database 
gets larger, the effects of decoherence on the efficiency of the quantum search algorithm cannot 
be ignored. Moreover, with decoherence, it may not be useful to iterate beyond the first maxima 
in the probability distribution of the search entry. Surprisingly, we also find that the number 
of iterations for maximum probability of the search entry reduces with decoherence. 

03.67.Lx, 89.70.+C 



I. INTRODUCTION 



In many computations, the search algorithm is one 
of the most time-consuming activity of the computer 
program |l|]. With the rapid proliferation of the in- 
ternet and the increasing need to retrieve information 
from this ever expanding internet network, it is crucial 
to consider alternatives to current search algorithms. 
Recently, quantum algorithms have been shown to re- 
duce hitherto computationally difficult problems for 
which there are no known polynomial time algorithm 
into a tractable problem involving polynomial time . 
This problem involves factoring a large composite in- 
teger into its prime factors. Besides prime number 
factorization, it has also been shown that quantum al- 
gorithms like the quantum search algorithm could sub- 
stantially improve the searching process for a database 
||-||] despite only a square root speed-up. 

Indeed, in recent years, this immense potential in 
quantum algorithm has spurred a fecundity of ideas 
on the actual implementation of a quantum computer. 
However, the prospect of actually building one in the 
next decade remains taunting and seemingly insur- 
mountable with current technology due primarily sev- 
eral reasons, like decoherence and scalability. Unlike 
classical bits, quantum bits (qubits) are highly sus- 
ceptible to collapse due to the difficulty of isolating 
the quantum mechanical systems from its environ- 
ment. Such decoherence inevitably leads to a loss 
of information within the system. Thus it is neces- 
sary to consider fault tolerant computation through 
quantum error correction [p| or more recently using 
decoherence- free subspaces |0, |j . 

Recently, some interesting works have been done to 
consider the robustness of the quantum search algo- 



rithm under a noisy environment ||,|T^. In ref. 
they modeled the decoherence of the search algorithm 
using a stochastic white noise. In this paper, by con- 
sidering a specific decoherence model, we obtain ana- 
lytic expression for the probability of the search item 
in a quantum search algorithm after a certain number 
of iterations. Our analytic expression can certainly 
facilitate the study of the behavior and efficiency of 
the search algorithm within the model. 

This paper is organized as follows. In section ||, we 
briefly describe the Grover's search algorithm under 
a noise- free environment. In section IE, we describe 



our decoherence model and in section [V derive an 
explicit analytical formula for the probability of the 
search item after a certain number of iterations within 
the model. In section |^, we consider the robustness 
of the search algorithm under the decoherence model. 
Finally in section ^l[ we summarize our results and 
discuss some implications of the model. 



II. GROVER'S SEARCH 



In a series of seminal papers, Grover |^-|5|] con- 
sidered a quantum algorithm that could achieve a 
speed-up in the computational implementation with 
only 0(a/]V) for a large structured database with N 
records. This algorithm compares favorably with the 
classical result which can only execute a search with 
0{N) efficiency. Moreover, it has been shown that 
Grover's algorithm is optimal 

Grover's search algorithm can be summarized 
neatly into the following main steps: (i) Initializa- 
tion of the system into a superposition of states; (ii) 
Subjecting the system to a hashing function, C(S), 



represented by a unitary operator, [/, given by 
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(assuming that the first entry satisfies the search cri- 
teria and therefore undergoes a rotation of tt radians) 
fohowed by a diffusion matrix, defined by 
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for 0{VN) iterations; (iii) Measuring the resulting 
state. The heart of the process therefore hinges sig- 
nificantly on the on step (ii) and the 0{^\/N) iterations 
of the matrix DU = S. In this paper, we look closely 
into this important step and scrutinized the ideas be- 
hind the efficiency. 

Before proceeding further, we first consider the 
eigenvalues of the matrix S. It is not difficult to show 
that the eigenvalues of S all lie on the locus \z\ = 1, 
unit circle, on the complex Argand plane and are ex- 
plicitly { — 1, • • • , — 1, 7^, 7^*} where the root r] and ry* 



(Ar-2) 

satisfy the equation z'^ 



2{N -2) 
TV 



z + 1 = 0. Indeed, 



it can be further shown that if the matrix S is diag- 
onalized as PAP~^, with A = Diag[-1, • • • , -1, r], r]*] 



(iV-2) 



and 



r]= — \N -2- 2iVN - 1 

then the matrix P assumes the form 
/ 



(3) 
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Since |?7| = 1, we can rewrite the complex number r] 
as e~^^, where, using eq(|l). 



N-2 



cost/ = 



N 



and 



(5) 



Finally, we note that the probability of finding the 
search item can be evaluated to be — (cos(m^) + 
\JN — \ sin(m^))^ while the probability of getting one 
of the other entries in the database is —{cos{mO) — 

, sin(m^))^. 



III. NOISE INDUCED THROUGH 
DECOHERENCE 



A. Superoperators 



Decoherence can be studied and understood within 
the context of superoperators through the evolution 
of a bipartite quantum system jl^]. A superoperator 
describes a linear map, $, from input density matrix, 
p^j^, to output density matrix, Pout- Although the su- 
peroperator, $, is not a unitary operator, it is a linear 
map that preserves hermiticity as well as the trace. 

The standard procedure of understanding the be- 
havior of one part of a bipartite quantum system is to 
extend the system to a larger one (in which the en- 
vironment (E) is incorporated) so that the evolution 
of state becomes unitary under the transformation, 
U. By assuming complete positivity of the superoper- 
ators, it is possible to study the non-unitarity evolu- 
tion of the state of a subsystem using an operator sum 
representation. In terms of the operator sum or Kraus 
representation, we can express this map, $, more suc- 
cinctly as 



(6) 



Unitarity of the operator, Uce^ also requires that the 
Kraus operators satisfy the condition 



J2m,mI = i 



(7) 



B. Decoherence Model 



For two- spin- 1/2 system, we can consider the fol- 
lowing Kraus representation 




where a\ i = x^y^z are the usual Pauli matrices. We 
can check easily that the unitarity condition for the 
larger space is satisfied since M^M^ = 1. This 



decoherence model is sometimes alluded to as the de- 
polarizing channel. 

This decoherence model is an idealization of a trans- 
mission and storage process with some elegant 
mathematical symmetry in which the quantum state 
in the channel can undergo a bit-flip and phase errors. 
The construction of a depolarizing channel arises from 
the interaction of a quantum state, |^) with an envi- 
ronment in which there is a probability p that the 
quantum state will survive and a probability p/3 that 
it would execute a pure bit-flip, a pure phase error or 
a combination of both. Thus, schematically, we have 





(12) 
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(15) 



Furthermore, it is possible to show that for general 
spin- 1/2 particle in which the density matrix describ- 
ing the mixed states is given by po = 1/2(1 + n • a), 
that the original density matrix po evolves under the 
depolarizing channel as 



1. 



(1 - P)PO 



(16) 



where I2 is the 2x2 unit matrix. In order to apply 
this model to the quantum search algorithm, it is nec- 
essary to consider the extension to N qubits. It is not 
hard to verify that in this case, Eq.([l6|) becomes 



p'=P-In + (1-p)po (17) 
where In is the N x N unit matrix. 



we allow the system to evolve under a depolarizing 
channel so that using Eq.([T7|), we have 



(20) 



It is easy to show inductively that for m iterations, 

Pm=P^lN(l + (l-p) + --- 

+{l-pr-^)+{l-prS'^poS^"'pi. (21) 



]^[i-(i-^'r] 



(1 - p)™5™po5'^" 



(22) 



A schematic representation of the model is shown in 
Fig. |l|. 
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FIG. 1. Schematic illustration of a decoherence model. 

After the m iterations, the state of the system 
is measured using the standard Von-Neumann mea- 
surement or positive operator valued measurement 
through appropriate extension to a higher dimensional 
space, so that the state after the measurement reads 



IV. SEARCH WITH DECOHERENCE 

Let us now study the effects of the decoherence (de- 
polarizing) model to the quantum search algorithm. 
As in any search algorithm, we first initialize the sys- 
tem into a superposition of states with equal proba- 
bilities. Thus the initial density matrix of the system 
is 



(18) 



This state is then subject to the usual inversion- 
diffusion process so that the density matrix at the end 
of the first transformation is 



Po 



SpoSl 



(19) 



In the noise-free Grover's search algorithm, this trans- 
formation is repeated a certain number of times. In 
our model, before we perform the second iteration. 



Pf = "^EiPmEj 



(23) 



where Ei are orthogonal projection opera- 
tors (POVM) However, we note that 
Y^^EiS'^poS'^^El is just the probability of finding 
the desired state under the search without decoher- 
ence. Moreover, using the results in section || or 
ref. P, 111, we have 



^EiS'^poS'^^E] = ^(cosmi9 + VaT^ sin m(9)^ 



(24) 



- . Finally, 



where cos6> = 1 — — and sin 6 

^ N N 
using Eq. ( |22| ) and manipulating using some simple al- 
gebra, we find that the probability of searching suc- 
cessfully for the desired state under decoherence is 



P(m) 



L{l-{l-pr[{cosme + 

VN - 1 sin m9 f 



1]) (25) 



Eq. ( p5D is the main result in our paper. 

For the case of a generahzed measurement, we 
can consider the set of r POVMs given by Fi = 
XijEj^ i = 1 • • • r, j = 1 ■ ■ ' N where the parame- 



ters Xij must satisfy the unitarity condition 



i=l 



1 for each j. For simphcity, we can chose r = N and 



(l-e), A, 



i j. Indeed, it is not 



difficult to imagine an experimental setup for these 
POVMs since they can be used to describe possible 
errors in the photon detectors. In this case, we find 
the probability of successfully searching for the desired 
state is 

P(m) = (l-e) P(m)ortho + 



— (1-(1- p)^ [(cos ml9 



: sin ml9)^ + 1]) (26) 



where P(^)ortho is the probability of the search item 
under Von Neumann orthogonal measurements and 
given in Eq. (p5|) . 



V. ANALYSIS 

Within the search algorithm, the number of itera- 
tions needed for the search is usually fixed by max- 
imizing the probability of getting the desired state. 
Without any decoherence, the number of iterations 
^max of the first maximum for large N is found to 
be 



^max 



(27) 



With decoherence, from Eq.(|25|), the condition is im- 
posed by the following transcendental equation. 



(l-y(l-cosC)) = ^smC 



1 — P 



(28) 



where C 2(2m + l)l9. 

It is interesting to compare the probability of the 
search item as a function of the number of iterations, 
m, for the specific case for various values of p. Such 
a graph for TV = 128 and p = 0.01,0.04 and 0.083394 
is shown in Fig. H. The value of p = Pc = 0.083394 
refers to the maximum allowed value of the decoher- 
ence parameter p subject to mmax = 7r\^/4 such 
that the probability of the search item do not fall be- 
low 0.5. Fig. H clearly shows that the presence of 



decoherence tends to provide a decaying effect on the 
periodic probability distribution of the search item. 
In the noise-free situation, we have iterations in which 
the probability of search items reaches a periodic max- 
imum without significant decay in the maximum prob- 
ability. In the case of a noisy search, we see that there 
is gradual reduction in the maxima so that subsequent 
maxima are usually rendered useless in the search as 
the magnitude of decoherence gets larger. 
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FIG. 2. Probability of search item as a function of the 
number of iterations for p = 0.1,0.4 and 0.83394 and 
N 128. 

As we have noted earlier, if we maintain the number 
of iterations to the value required under the noise-free 
situation, namely mmax, we can solve Eq.(p5|) for the 
value of Pc such that the probability of the search item 
is above chance level (i.e. P(m) > 0.5) for p < Pc. In 
Fig ^, we plot the critical decoherence parameter p as 
a function of the database size n = log2{N) up to a 
database size of about 2 x 10^. This graph shows that 
as the database gets larger, it is imperative and nec- 
essary to overcome the problem of decoherence. At a 
database size of 2^^ ^ 2x 10^, the amount of decoher- 
ence allowed is only about 0.000609, which is a very 
small value indeed. 
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FIG. 3. Critical decoherence parameter pc as a function 
of n = log2{N) 

In a general search problem, the number of iter- 
ations can always be controlled. For a noisy search 
algorithm, it seems from our analysis that it may be 
better to reduce the number of iterations so that the 



probability of the search items becomes a httle larger. 
Using a database of = 1024, we plot the variation 
of the probability of the search items as a function of 
iterations m for p = 0.001, 0.014 and p = Pc = 0.0274 
in Fig. ^. The graph clearly shows a shift in the max- 
ima for the curves. Indeed for the near critical value 
of p = Pc, we see that the maxima hovers around the 
value of 0.5 at m = 21. A three-dimensional plot of 
the probability as a function of m and p is shown in 
Fig ||. Indeed, by varying the decoherence parameter 
p, it is possible to solve for the value of mmax corre- 
sponding to the probability of the search item. Fig. || 
shows the value of mmax as a function of p. The step 
function refers to the integer part of mmax- 
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FIG. 4. Probability of the search item as a function 
of the number of iterations m for p — 0.001,0.014 and 
p = Pc = 0.0274 






FIG. 5. Three-dimensional plot of the probability of the 
search item as a function of the number of iterations, m 
and decoherence parameter, p (N = 1024) 




0.005 0.01 0.015 0.02 0.025 

P 

FIG. 6. The value of mmax corresponding to the max- 
imum probability of the search item as a function of p 



It is interesting to note that the graph of mmax cor- 
responding to the maximum probability of the search 
item is a linear function of p. A interpolation fit gives 
mmax = 24.6254 - 127.7426p. If we then plot the 
probability of the search item as a function of p using 
the integer part of this linear function, z.e.Int(mmax) 
and compare it with the corresponding probability for 
^max = Int(7rV^/4), we see that the difference is 
reasonably negligible, especially for small p. This plot 
is illustrated in Fig. |^. 
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FIG. 7. Probability of the search item as a function of 
p for different iteration formulae. 



In the context of a generalized measurement, we 
find that probability amplitude of the search item is 
given by the expression in Eq. ( p6|) and this probability 
generally reduces in the first maximal iteration. The 
plots of the probability of the search item as a func- 
tion of the number of iterations m for various values 
of j9 at = 128 and r = 131 is plotted in Fig. |. 
In the plots, we have chosen a reasonably small value 
for the parameter, e = 0.1, since under experimental 
conditions, one do not expect a large value for this 
parameter. In fact, the results show that for suffi- 
ciently small e, it is possible to confine our study to 
orthogonal Von Neumann measurements. 
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FIG. 8. Probability of the search item as a function of 
m for iV = 128, e 0.1 and (a) p 0.01, (b) p 0.04 and 
(c) p = 0.083394 for orthogonal (bold) and non-orthogonal 
(dotted) measurements. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 



In this paper, we study the effects of the search al- 
gorithm analytically and numerically using a simple 
decoherence model. Several observations are in order. 
Firstly, we see that the search algorithm can be ren- 
dered useless for small degree of decoherence under 
the depolarizing channel as the database gets larger. 
Indeed for a database size of about 2 million entries, 
the degree of decoherence, p, is a factor of 1000 smaller 
than the corresponding value for a database size of 4. 
It appears that for large A/", the factor varies approx- 
imately as a/]V. 

Secondly, as we would expected, decoherence gen- 
erally diminishes the probability of the search items. 
With high decoherence, it may not be useful to iterate 
beyond the first maxima in the probability distribu- 



tion of the search item. 

Thirdly, the number of iterations needed to maxi- 
mize a search probability decreases with increasing de- 
coherence. This is a somewhat unexpected result since 
decoherence ostensibly tends to reduces the probabil- 
ity distribution. Nevertheless, although the associated 
maximum probability has decreased, it is interesting 
to note that under a noisy channel, one could in gen- 
eral improve the efficiency by reducing the iterations 
prior to measurement. It is interesting to note here 
that the number of iterations (considered as a contin- 
uous variable, i.e. prior to taking the integer part) 
falls off linearly with the degree of decoherence, p. 
However, if we continue to use the maximum number 
of iterations in noise-free situation, the difference in 
probability appears to be less than 10%. 

Finally, we note that although the decoherence 
model that we have employed here is very simple, it 
certainly provides a wealth of information regarding 
the behavior of the search algorithm under a noisy 
environment. 
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